Reply to "Incommensurate vortices and phase transitions in two-dimensional XY 
models with interaction having auxiliary minima" by S. E. Korshunov 



o 

(N 

^1" 



o 
o 



cn 
o 



X 




Gabriel A. Canova, 1 Fabio Poderoso, 1 Jeferson J. Arenzon, 1,2 and Yan Levin 1 

1 Instituto de Fisica, Universidade Federal do Rio Grande do Sul CP 15051, 91501-970 Porto Alegre RS, Brazil 
2 Universite Pierre et Marie Curie-Paris VI, LPTHE UMR 7589, 
4 Place Jussieu, FR-75252 Paris Cedex 05, France 

We present a rigorous proof and extensive numerical simulations showing the existence of a tran- 
sition between the paramagnetic and nematic phases, in a class of generalized XY models. This 
confirms the topology of the phase diagram calculated by Poderoso et al. [PRL 106(201f )067202]. 
The results disprove the heuristic argument presented by Korshunov in arXiv:1207.2349vl, against 
the existence of the generalized-nematic phase in a model with q — 3. 



In a recent Letter [![, we have studied the phase 
diagram of a generalized XY model with Hamiltonian 
H = -£ fo )[Acos(0 4 - %) + (1 - A)cos(<^ - #,)], 
where q > 1 is an integer and < A < 1. Using Monte 
Carlo (MC) simulations, we showed that for q = 3 and 
< A < 0.4, the model exhibits — depending on the 
temperature — three possible phases: paramagnetic (P), 
generalized-nematic (N), and ferromagnetic (UF). The 
phase transition between P and N was found to be in 
the Kosterlitz-Thouless (KT) universality class, while the 
transition between N and UF was found to belong to the 
3-state Potts universality class. In his Comment on our 
work [2]. Dr. Korshunov argued that the N phase cannot 
exist for q > 2, and that there should only be one "gen- 
uine phase transition" between the P and UF phases. We 
will now show that the argument of Ref. [2j is incorrect. 

Let us first consider A = 0. In this case the Hamil- 
tonian becomes purely g-nematic. Changing variables in 
the partition function, qOi — > 6i, shows that the model 
is isomorphic to the usual XY model, but with the low 
temperature phase N, instead of UF. The phase transi- 
tion from P to N will, therefore, occur at To ~ 0.893, 
the same temperature as for the standard XY model and 
will belong to the KT universality class. This, clearly 
demonstrates that the N phase exists for A = 0. Using 
Ginibre's inequality [|[ it is possible to show that the P to 
N transition will also extend to finite A Q • Furthermore, 
Ginibre's inequality allows one to derive a rigorous lower 
bound |4J on the transition temperature between P and 
N phases, Tkt(A) > (1 — A)T . Since at very low tem- 
perature the system must be in UF phase, this proves 
the existence of P, N and UF phases for small, but fi- 
nite values of A. contradicting the heuristic argument of 
Ref. 0]. 

To precisely delimit the location of all three phases for 
the model with q = 3, we consider a specific example, 
A = 1/4. For this A, and using finite size scaling (FSS), 
in Ref. [l| we have calculated the critical temperature 
for the N-UF transition to be Tp ot ts — 0.365, which was 
found to belong to the 3-state Potts universality class. 
The order parameter mi (magnetization) shows clearly 
this transition, see Fig. Q] On the other hand, at Tkt, 
the nematic order parameter 7713 shows the transition 
between N and P phases. At the transition tempera- 



FIG. 1: Order parameters mi and ms (see Ref. PI) for several 
system sizes L showing phase transitions at Tp tts — 0.365 
and Tkt — 0.68. Inset: Helicity modulus T versus T. The 
crossing with the line 2T/ir at Tkt(£), when extrapolated to 
L — > 00, gives Tkt — 0.68. 



ture, 7773 decreases with L as 7?73(Tkt) ~ L^^l v ', with 
(3/v ~ 0.117. The exponent is very close to the theoret- 
ical value expected for the KT transition, 1/8 = 0.125. 
To further verify the "genuineness" of this transition, we 
calculated the helicity modulus T Q , shown in the inset 
of Fig. [T] as a function of temperature, for several sys- 
tem sizes. The helicity modulus crosses the straight line 
2T/-7T [6] at Tkt(L) and, extrapolating to L — > 00, we ob- 
tain Tkt — 0.68, slightly above the lower bound provided 
by Ginibre's inequality. 

In Fig. [2] we present the susceptibility \3 as a function 
of T for different system sizes. The phase transition is 
very clear from the divergence of the susceptibility at 
Tkt, as L — > 00. For a KT phase transition, the FSS 
predicts that X3(^kt) ~ i 1 ' 75 , while our simulations find 
£-1.766 Fina n yi if we p i ot X3 L 7 '~ 2 , with the KT 77 = 

1/4, vs. the Binder cumulant, all the susceptibilities for 
different system sizes should collapse onto a universal 
curve 8]. This is precisely what is found in our MC 
simulations, see inset of Fig. [2J 

Ref. 0] also questions the transition between the 
phases Fi and UF, in the model with q = 8, and the ab- 
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sence of any "qualitative" difference between these two 
phases. This, however, is clearly not an issue, as is ex- 
emplified by the usual liquid-gas phase transition — the 
difference between liquid and gas being only "quantita- 
tive" . The Fig. 6 of Rcf. [1] shows clearly the transition 
between Fi and UF, which again belongs to the KT uni- 
versality class. 



T 



FIG. 2: The susceptibility \3 associated with mz near the KT 
transition for various system sizes. Inset: rescaled suscepti- 
bility versus the Binder cumulant, Us = (m^) 2 / (m|) , showing 
a perfect collapse with the KT exponent r\ = 1/4. A similar 
collapse is also obtained for the magnetization 



In conclusion, we have presented a rigorous proof, as 
well as, numerical evidence for the existence of a transi- 
tion between the N and P phases belonging to the KT 
universality class, at odds with the heuristic argument of 
Rcf. 0. 
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